Local Gaussian operations can enhance continuous- variable entanglement distillation 
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Entanglement distillation is a fundamental building block in long-distance quantum communi- 
cation. Though known to be useless on their own for distilling Gaussian entangled states, local 
Gaussian operations may still help to improve non-Gaussian entanglement distillation schemes. 
Here we show that by applying local squeezing operations, both the performance and the efficiency 
of existing distillation protocols can be enhanced. We derive the optimal enhancement through local 
Gaussian unitaries, which can be obtained even in the most natural scenario when Gaussian mixed 
entangled states are shared after their distribution through a lossy-fiber communication channel. 
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Entangled quantum states can be seen as one of the 
essential resources in quantum information processing 
(QIP) [1]. However, entanglement is fragile and easily 
degraded by uncontrollable environment-induced noise, 
for instance, during its distribution over longer distances. 
In order to circumvent this problem, various entangle- 
ment distillation protocols have been proposed capable of 
improving shared entanglement by local operations and 
classical communication [2] . In the realm of continuous- 
variable QIP, although Gaussian entangled states and 
Gaussian operations [3, 4] are widely used for telepor- 
tation, dense coding, and other protocols [5], distilling 
Gaussian entanglement in the potentially most efficient 
way through only Gaussian operations has been shown 
to be impossible [6-8]. 

In the mean time, efforts have been made to incorpo- 
rate the necessary non-Gaussian element into optical en- 
tanglement distillation schemes. For example, Opatrny 
et al. introduced the photon subtraction (PS) strategy 
employing one of the most readily available non-Gaussian 
operations to distill stronger entanglement [9]. This 
type of distillation was recently demonstrated experimen- 
tally [10]. In other, complementary experiments, non- 
Gaussian entangled states were initially prepared and 
subsequently, since no longer restricted by the above no- 
go results, distilled through Gaussian operations [11, 12]. 
All these results and developments seem to imply that 
local Gaussian operations are useless for the distillation 
of Gaussian entangled states (see ,for example, Ref. [13]; 
however, note [14]). 

In the present work, we shall demonstrate that Gaus- 
sian operations, though insufficient for distilling Gaus- 
sian entanglement on their own, can still be used to en- 
hance existing non-Gaussian distillation schemes. More 
specifically, we present examples of a distillation protocol 
that makes use of local squeezing operations in addition 
to Opatrny et al. 's non-Gaussian PS strategy. We then 
show that after distillation, both the success probability 
and the entanglement (in terms of logarithmic negativity) 
are improved in a regime when an initial Gaussian two- 



mode squeezed vacuum state (TMSS) is weakly squeezed. 
Moreover, this result even holds when the TMSS is first 
subject to an amplitude damping channel, as it would 
be the case, for instance, during an optical-fiber-based 
entanglement distribution. In fact, we are able to show 
that the local Gaussian pre-processing prior to PS-based 
distillation (but after the channel transmission) may be- 
come even more important for a lossy channel. 

Our result is of both conceptual and practical signifi- 
cance. From a more fundamental point of view, it shows 
that in order to optimize quantum information tasks as 
important as entanglement distillation, discrete-variable 
and continuous-variable techniques should go hand in 
hand [15]. Since the PS technique [16] as well as the im- 
plementation of online squeezing [17] is becoming state 
of the art, our protocol is of practical relevance too. 

Preliminaries- for representing a continuous-variable 
system consisting of N bosonic modes, we shall employ 
the phase-space description, where every mode k can 
be conveniently expressed by field quadrature position 
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FIG. 1. Distillation of Gaussian continuous-variable entan- 
glement through PS and local squeezing, Sa(t~'a) an d Sb(t'b)- 
Here, \iP)ab is a pure TMSS. Beam splitters with transmis- 
sion coefficient T and on-off detectors C and D are used to 
realize the corresponding non-Gaussian PS operations. 
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and momentum operators, 
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(fife — fi fe )/(i\/2), with fife, fi fe being the mode annihi- 
lation and creators operators. Throughout, we use 
hats to denote operators in Hilbert space. The canon- 
ical commutation relations for an TV-mode system can 
be conveniently written as [X m ,A"„] = ifi„ m , 1 < 
m,n < 2N, with X = (xi,p\, ■ ■ ■ , Xn,Pn) and = 
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Furthermore, the density matrix of 
p in infinite-dimensional 
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the above iV-mode system 
Hilbert space is represented by the characteristic func- 
tion x(0 = Tr{pexp[?X T ^]}, £ € R 2N , or equivalently by 
its Fourier transform, i.e., the Wigner function W(x) = 



■ exp [— ix T £] x{0- I R particular, Gaussian 
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states are those states whose characteristic or Wigner 
functions are Gaussian: x{0 — ex P [ — 5£ T ^£ + * xT £] , 

and W(x,x,V) = cxp [-y x ~y^ (x - g) ] . Such Gaus- 

sian states are fully characterized by displacements x = 
Tr(Xp) and a covariance matrix (CM) V, with entries 
Vi m = \{XiX m +X m Xi) - (Xi)(X m ). For example, for 
the pure TMSS, we have x = and the CM is given by 



V AB = - 
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with c = cosh2r, s = sinh2r. Unitary state evolutions in 
Hilbert space generated by Hamiltonians quadratic in the 
canonical operators correspond to symplectic transforma- 
tions, S € Sp(2N,R), in phase space. In what follows, 
we will frequently use the symplectic single-mode squeez- 
ing, S(r) = diag{e r , e _r }, and two-mode beam-splitting 
operations. The latter is performed via the 4x4 matrix 
Sgg with non-vanishing elements (5*bs)" = Vi, and 
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Here, T represents the transmission coefficient and kl the 
two input modes of the beam splitter. 

Pure-state distillation- with the notations above, we 
can now proceed to derive the state evolution in our 
modified entanglement distillation protocol, as shown in 
Fig. 1. As a start, let us first assume that the initial state 

oo 

is the pure TMSS, \ip) AB = £ Vl - X 2 X n \n) A \n) B , A = 

71=0 

tanhr, with \n) denoting the Fock basis and r the squeez- 
ing parameter. To perform entanglement distillation, we 
propose to use two local squeezing operations S A (r' A ) and 
S B {r' B ) (with squeezing parameters r' A and r' B ) before 
the PS operation. For the latter, we employ a beam split- 
ter (with transmission coefficient T) and conventional on- 
off detectors [18]. Throughout this paper, we refer to a 
successful distillation event when both detectors register 
nonzero photon counts. 

Now assume modes C and D are initially prepared in 
a vacuum state. Thus, the CM of the four-mode state 



ABCD (before the photon detections) becomes 
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[S A (r')®S B (r')], (3) 



under the additional assumption that optimal distillation 
will occur for an initial symmetric state through symmet- 
ric local PS and squeezing operations, r' A = r' B = r' . 

In our scheme, we propose to use on-off photon de- 
tectors, as commonly employed in quantum optics ex- 
periments. Such a detector is represented by two mea- 
surement outcomes: 'off' when no photons are de- 
tected and 'on' when one or more photons arc de- 
tected. Through a successful distillation event, modes 
C and D are projected onto non-vacuum components 
and the state of modes A and B is reduced to p = 



Tr C£ , 



o i <x> n(°™) <x> ri(°™) 

r ABCD AB ^ ±± c W D 



j /p succ , with n(°") = 

0) (0 1 = Y^=i \ n )( n \- Throughout, we use I m to 
represent an m-dimcnsional identity matrix. In order 
to obtain analytical results, we shall again employ the 
phase-space formalism. In fact, although the single-mode 
operator n(°") leads to a non-Gaussian Wigner function, 
W(x) = £~ ~ \ cx p[~ x T /x] [19], by expressing every 
single operator through a Wigner function and carrying 
out the corresponding integrals, we find that the Wigner 
function of the distilled state is a linear combination of 
four Gaussian functions: 
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with Pi = 1, P 2 = -[dct(V c + / 2 /2)]- 1 / 2 , P 3 = 
-[det(V D +h/2)]- 1 '\ Pi = [det(r C D+V2)]- 1/2 , 
Vi = T AB , V 2 = T AB — <?\ (V c + h/^)' 1 aJ,V 3 = T AB - 
02 (V D + h/2)- 1 o-J, V 4 = Tab o~ (T C d + h/2)- 1 a T , 
and Tab,Tcd, Vc, Vd defined by partitioning Vabcd as 
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Here, a = ( o\, d 2 ), and o\ and a 2 are both 4x2 ma- 
trices, T A b,Tcd are 4x4 matrices, and Vc,Vd,C are 
2x2 matrices. By also taking into account the nor- 
malization of the Wigner function and integrating both 
sides of Eq. (4) over the whole phase space, the proba- 
bility of successful distillation can be found to be given 

by P succ = Ej=i p j- 

Let us now use the logarithmic negativity [20-22] as a 
figure of merit to quantify the entanglement of the dis- 
tilled state. The logarithmic negativity of a bipartite 
state p ab i s defined as Ptv(pab) = log 2 HpabII w ^h Pab 
being the partially transposed density operator. In order 
to compute the logarithmic negativity, we shall calculate 
the density matrix from notions in phase space. Using a 
method similar to Ref. [23], we have the following theo- 
rem [24]. 
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FIG. 2. Distillation of a pure TMSS (r = 0.025) using local 
squeezers (r') and PS method (T = 0.95). (a) logarithmic 
negativity of the output state; the blue circle indicates the 
optimal r' which maximizes the logarithmic negativity (ropt = 

0. 1565, En = 1.5278) (b) success probability of distillation, 

1. e., the probability that both detectors obtain 'on' results. 



Theorem 1: let V be the CM of a two-mode zero- 
displacement Gaussian state and pab be the correspond- 
ing density matrix in the Fock basis, then the normalized 
matrix elements (for any k\,k2,m\,m,2 — 0,1,2,...) fol- 
low from 
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where o~ x ,cr z are the usual 2x2 Pauli matrices. 

For our non-Gaussian state whose Wigner function 
is a linear combination of Gaussian functions, VF(x) = 
J2j :PjW(x, 0, Vj), it can be easily proved that the den- 
sity matrix follows the linear rule p = ^ . Pjp(Vj), with 
p{Vj) determined by Theorem 1. 

Using the above methods, we now present a numerical 
evaluation of the performance of our modified distilla- 
tion protocol. Fig. 2 shows the logarithmic negativity 
of the distilled state for r = 0.025, r' £ [0.01,0.20] [25]. 
The logarithmic negativity attains a maximal value at 



an intermediate point r' opt = 0.1565, whereas the success 
probability increases monotonically with local squeezing 
r' . At the optimal point r' t , with success probability 
3.5029 x 10" 6 and E N = 1.5278, we obtain a signifi- 
cant improvement over Opatrny 's original PS strategy 
(with 1.5645 x 10~ 6 and E N = 0.1352 [26, 27]). In other 
words, the local filter operations corresponding to the 
beam splitters and on-off detectors (i.e., the local map 
for say mode A, £ : p A -> J^Zi ^PaEJ /Tr^ e\e iPa ], 

with Ei = £^i(-i)*>/(5T(*»-*>/2(i - r) l /2|„ _ i)( n \) t 

become more efficient when replaced by the set of opera- 
tors {EiSA(r' A )}i=i,2---- At the same time, local squeez- 
ing increases each mode's average photon number and 
thus enhances the probability that photons are detected. 

Mixed-state distillation- in the above distillation 
scheme, the local squeezers could be as well seen as part 
of the initial state preparation. One could then argue 
that so far we have simply found the optimal Gaus- 
sian pure state for PS-based entanglement concentration 
which turns out to be different from the TMSS. It is 
therefore intriguing to examine whether the TMSS would 
also benefit from local Gaussian pre-processing after its 
transmission through an imperfect channel such as a lossy 
fiber. This would lead to a clear distinction between lo- 
cal operations before and after the entanglement distri- 
bution, corresponding to alternate state preparations or 
alternate state distillations, respectively. 

For this purpose, we use beam splitters and auxiliary 
vacuum modes to model the optical amplitude-damping 
channel [27]. Our method introduced above still applies, 
except for replacing c — > c' = 1 — r/ + r) cosh 2r, s — ¥ 
s' = r)smh2r in Eq. (1), with rj being the transmis- 
sion efficiency of the lossy channel. We find that local 




FIG. 3. Mixed-state distillation with 77 = 0.5. Comparing 
the distillation performances between Opatrny's PS strategy 
(blue, dashed line) and our local-squeezing enhanced PS strat- 
egy (red, dot dashed line), for r G [0.005,0.4], T = 0.95, 
r' — r = arctanh(A). The black (solid) line indicates the 
logarithmic negativity of the mixed state before distillation. 
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FIG. 4. Teleportation fidelity, (a) with n = l,r = 0.025 (b) 
optimal local squeezing r' optF as a function of r (c) comparison 
between our modified scheme (red, dot dashed line) and Opa- 
trny's PS strategy (blue, dashed line): the results coincide 
beyond r ~ 0.13 (d) using a distilled 3dB-amplitude-damped 
TMSS state; the threshold value becomes r ~ 0.26. The black 
(solid) line shows the fidelity of using the 3dB-damped TMSS 
state before distillation. Throughout we use T = 0.95. 



squeezing still helps to improve entanglement distillation 
of mixed states. In Fig. 3, the logarithmic negativity 
and success probability of distillation are shown. For 
simplicity, we consider the case when each mode of the 
TMSS is transmitted through a 3dB (77 = 0.5) amplitude- 
damping channel. In the low-squeezing regime, e.g. for 
r E [0.005,0.4], a significant improvement is obtained. 

Quantum teleportation- the improvement in our mod- 
ified scheme can be experimentally verified in an oper- 
ational fashion through quantum teleportation [9, 26]. 
Let us consider standard unit-gain teleportation in which 
the entangled state after distillation is pre-shared and 
an unknown coherent state |ao) is to be teleported. 
The teleportation fidelity is _F t ci = (ao|y out(Q!o)|Q ! o), 
with Pout(cKo) being the ensemble average of all output 
states conditioned upon different Bell measurement re- 
sults. Due to its linearity, the teleportation fidelity us- 
ing the entangled state in Eq. (4) can be conveniently 



written as F = Y,j=i P jFtei(Vj)/P SUC c, with (for any 
ao) *tei00 = det(RaR + 7 T i? + Rj + + I 2 y 1/2 ,R = 

( Ot / Y 

diag(— 1, 1), where a, /?, 7 are defined by V = ( ^ T ^ 

In Fig. 4(a), we consider a pure TMSS state (77 = l,r = 
0.025) being distilled and then used for teleporting an un- 
known coherent state. Fig. 4(b) shows the optimal local 
squeezing r' optF as a function of r. Clearly, r' optF begins 
to drop towards zero for a stronger pure TMSS state 
(r ~ 0.13), representing a threshold beyond which our 
modified scheme ceases to improve Opatrny's PS strat- 



For a 



shows the fidelity for using a distilled 3dB-damped TMSS 
state. In this larger threshold value, r ~ 0.26, 

occurs. Thus, we find that for an amplitude-damped re- 
source state, the local Gaussian pre-processing prior to 
distillation becomes even more significant. 

Optimal Gaussian unitaries- the most general sym- 
plectic transformation applicable to our local filters 
can be written as a sequence of phase rotation, local 
squeezing, and another phase rotation [28], U(r,0,<f>) — 
cos 8 sin 8 \ / e r \ / cos cf> sine; 
-sin0 cos 8 J \ e~ r J ^ — svaip cos< 
fixed squeezing r = arctanh(A) in the initial pure TMSS 
of Eq. (1), we applied two local unitary operations 
U (r, 8a , 4>a ) and U(r,0B,<l>B) to modes A and B , respec- 
tively. We computed the logarithmic negativity and suc- 
cess probability for 3000 randomly chosen 8 'a,4>a,8b,4>b, 
and we found that both figures of merit attain maximal 
values for 8a = 4>a = 83 — §b = 0. 

Summary- we demonstrated that local Gaussian oper- 
ations can be useful to enhance the distillation capabil- 
ities of non-Gaussian operations applied upon Gaussian 
entangled states, even when the initial states are Gaus- 
sian mixed states such as those emerging from an imper- 
fect channel transmission for realistic entanglement dis- 
tribution in quantum communication. In our protocol, 
for the distribution of one entangled-state copy, we con- 
sidered a photon loss channel, and for one-copy distilla- 
tion, we used the optimal local unitary squeezers in addi- 
tion to photon subtraction. It remains an open question 
whether local Gaussian non-unitary maps could further 
improve non-Gaussian entanglement distillation schemes. 
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e gy i r ' = 0)- I n Fig- 4(c), the fidelity at r' 
is compared with Opatrny's PS strategy, while Fig. 4(d) 
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